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Abstract 

The standard kicked top involves a periodically kicked angular momentum. By considering this 
angular momentum as a collection of entangled spins, we compute the bipartite entanglement 
dynamics as a function of the dynamics of the classical counterpart. Our numerical results indicate 
that the entanglement of the quantum top depends on the specific details of the dynamics of the 
classical top rather than depending universally on the global properties of the classical regime. 
These results are grounded on linking the entanglement rate to averages involving the classical 
angular momentum, thereby explaining why regular dynamics can entangle as efficiently as the 
classically chaotic regime. The findings are in line with previous results obtained with a 2-particle 
top model, and we show here that the standard kicked top can be obtained as a limiting case of 
the 2-particle top. 

PACS numbers: 03.67.Bg,05.45.Mt,03.65.Sq 
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I. INTRODUCTION 



The quantum-classical correspondence is the hallmark of semiclassical systems. These are 
genuine quantum systems for which the semiclassical expansion in the path integral prop- 
agator holds (at least for some appropriately chosen dynamical or temporal regimes). The 



le properties of a quan- 
ll. This is particularly 



quantum-classical correspondence allows to compute and interpret t 
tum system in terms of the properties of its classical counterpart 
important for systems displaying a complex dynamics, for which exact quantum computa- 
tions are either unfeasible or yield numerical results that hardly give any clues allowing to 
grasp the dynamics of the quantum system. 

These last few years several studies aiming to apply the quantum-classical correspondence 
to the understanding of dynamical entanglement have been published 2|-l8|]. Although en- 
tanglement is a distinctive quantum feature without a classical counterpart, many quantum 
systems displaying entanglement have a classical counterpart. The idea in these type of 
studies is then to assess whether there is a link between the generation of entanglement in 
the quantum system and the underlying classical dynamics. It was initially suggested that 
underlying chaotic dynamics was correlated with higher and faster entanglement. However 
it was later realized that integrable dynamics could lead to equivalent or even more efficient 
entanglement 

In earlier works , we have shown that the entanglement dynamics did depend on 
classical phase-space features, but in a specific and system-dependent way rather than in 
a generic manner. Our results were obtained in a modified two-particle kicked top, i.e. 
a kicked top involving explicitly two coupled angular momenta, one for each of the two 
entangled particles. This modified kicked top is different from the standard kicked top, 

n 

which is well-known to be one of the prototypical systems of quantum chaos [10[. Indeed, 
the standard kicked top involves a single angular momentum kicked by an external force 
whereas in the model we employed, an angular momentum is kicked by an interaction with 
the second particle, inducing a change in both particles' angular momentum. Accordingly, 
the classical dynamics of our modified kicked top is richer than the dynamics of the standard 
top. Notwithstanding, it can be rigorously shown (see below) that the standard kicked top 
is a limiting case of the modified kicked top we employed in our previous studies. 

In the present paper, we will study the entanglement generation as a function of the un- 



derlying classical dynamics of the standard kicked top. Although the standard top involves a 
single angular momentum, this angular momentum can be considered as a composite angu- 
lar momentum resulting from the entanglement of several elementary angular momenta (eg, 
qubits). Previous results lll-ll3| on the standard top employed in this context suggested 



that entanglement generation in the quantum top is correlated with chaos in the classical 
counterpart, and that in this system entanglement could be seen as a signature of chaos. 
These results appear to conflict the results we obtained with the modified kicked top. The 
physical origin of entanglement is surely different as one goes from the modified to the stan- 
dard top, but the nature of the quantum-classical correspondence is not expected to change. 
We therefore investigate in this work the relationship between entanglement generation and 
the quantum-classical correspondence in order to assess to what extent entanglement can 
be taken ClS db SI gnature of chaos. 

We will start in Sec. II by introducing the standard kicked top, not in the usual way 
but as the limiting case of our modified kicked top employed in earlier works. By doing 
so we will establish the relationship between these two models both based on stroboscopic 
maps. We will then describe the single standard kicked top as a compound system (Sec. 
Ill); we will see that the reduced linear entropy (that we will take to be the marker of the 
entanglement rate) depends on the averages of the angular momentum projections. The 
entanglement rate is maximized when the sum of these averages is minimized. From the 
quantum-classical correspondence viewpoint, the pertinent variable will consist in obtaining 
the classical dynamics minimizing this sum. The results are given in Sec. IV, for regular, 
mixed phase-space and chaotic dynamics. We will see that generically, chaos indeed mini- 
mizes this sum, but for appropriately chosen initial states, regular dynamics can entangle 
more efficiently, and in a more controlled fashion. We will discuss our results and conclude 
in Sec. V. 



II. THE KICKED TOP: AN ALTERNATIVE DERIVATION 



We present in this section t 



in our previous studies 



0, I9 



le link between the Rydberg molecule model, that we employed 



14| on entanglement generation and the quantum-classical 



correspondence, and the well-known kicked top (whose entanglement properties relative to 
the underlying classical dynamics will be studied in Sees. III-V; readers solely interested in 



the kicked top results may jump directly to Sec. III). The aim of this section is to show that 
the standard kicked top can be seen as the limiting case of the Rydberg molecule model 
when the total angular momentum and one of the two coupled angular momenta become 
infinitely large. 



A. The Rydberg molecule model: torsion and rotation of two coupled angular 
momenta 

1. Historical Introduction 

Electronic states of atoms or molecules are called Rydberg states, as opposed to valence 
states, when an outer electron moves far away from the remaining ionic core. These states 
form electronic series which converge towards the ionization limit of this outer electron. 

The starting point of the quantum analysis of such states was the Quantum Defect Theory 
(see e.g. the review article by Seaton [la]), established first for atoms. It was shown that, 
due to the non zero spatial extension of the ionic core, the levels near the ionization limit 
follow the hydrogen Rydberg law En = — Ry/(n + d)'^, with only a constant (or nearly so) 
shift d of the principal quantum number n, entitled Quantum Defect. Quantum Defect 
Theory was extended to Multichannel Quantum Defect Theory (MQDT), for the case that 
there are several series which converge to nearby states of the ion, and interact strongly. 
This theory depends only on a small number of parameters, basically one quantum defect 
per interacting series. Practically, all is solved with matrices whose size is the number of 
series, while "brute force" methods would in principle try to diagonalize a matrix which 
contains an infinite number of levels for each series. 

This theory was extended to molecules by Fano 
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171]. There are always many in- 



teracting series corresponding to the rotational states of the ionic core. Indeed the slow 
velocity of the core rotation leads to a splitting of the rotational states of the core which 
is of the same order of magnitude as the splitting between high lying electronic Rydberg 
states. The novelty was the implications of the anisotropy of the core. The effect of this 
anisotropy on the ionic potential decays faster with distance r than the point charge 1/r 
Coulomb potential, at least as 1/r^ or 1/r^. Fano showed that the key point of the analysis 
is the existence of a cut off distance Tq. Below this distance the motion of the outer electron 



is tightly bound to the direction of the ionic core, above it the two become independent. 
Many detailed studies have followed on moderately excited Rydberg states of molecules, see 
e.g. reviews in refs. [l8|, Il9 |. 



2. Phase Space: Dimension and Coordinates 



For a diatomic molecule this problem is in principle a three body problem, the two ions 
which constitute the molecular core and the Rydberg electron. After separation of the center 
of mass motion it depends in configuration space on six parameters, which can be chosen as 
the coordinates R,9mi^m of the relative position M of the two ions, and the coordinates 
r, 6'e, v^e of the Rydberg electron (all in the laboratory frame moving with the center of mass). 
The two coupled angular momenta we will study here are the angular momentum of the core, 
N, associated with the angles 6m, ^m-, and the angular momentum of the Rydberg electron 
L, associated with the angles d^.^e (see Fig. [1]). 

The problem for our purpose is simplified by two approximations 
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2l(, and the 



modulus L of L are kept constants. The first amounts to neglecting the vibrational motion 
of the core, and is frequently fairly well valid, since vibrational quanta are much greater 
than rotational quanta of the core. The second supposes that the isotropic part of the 
non Coulomb short range part of the Rydberg electron - ionic core interaction potential is 
much greater than its non isotropic part. Its validity is often acceptable. With these two 
approximations, the dimension of the classical phase space is decreased from 2 * 6 by 2 * 2 to 
8. Taking into account the conservation of the total angular momentum J and its projection 
Jz on the laboratory axis, decreases further this dimension by 2 * 2, i.e. the dimension is 
equal to four. Two classical phase space coordinates can be chosen as r and its conjugate 
momentum p^- For reasons which come from the quantum treatment below, the other two 
can be chosen either in the collision frame as the projection Lzq of L onto the OZq = M 
axis and its conjugate angle, the angle of the projection of L into the OXqYq plane with 



the OXq axis (known as 7r/2 plus the ascending node in classical mechanics treatises 22|). 



or in the laboratory frame as the modulus of the momentum of the ionic core and its 



associated angle ip^, which precise value is given in [2l|, Eq. (A. 11)] 




FIG. 1. (a): The laboratory frame, in which the core (depicted by the the two atomic nuclei) 
rotates approximately freely. M is along OZq. (b) The collision (or molecular) frame: it is the 
preceding OXqYqZq rotated around OZq so that the new OX is along N. During the collision 
L rotates around OZq by an angle 5ipL. The "free" rotation appears in this frame as a rotation of 
L around the new OX. 

3. Quantum solution: MQDT as a quantum map 

The usual manner of obtaining the solutions with MQDT is described in the Appendix. 
It involves introducing two different angular bases, the collision basis and the free rotation 
basis. The former gives the correct physical description when the electron is near the molec- 
ular core (and therefore collides) whereas the latter corresponds to a large radial distance 
between the electron and the core (which then rotates freely in the laboratory frame). For 
our present purposes, it is however more meaningful to envision the quantum problem as a 
quantum map. Indeed the quantization condition given by Eq. (lA.lOp obtained from MQDT 



can be rewritten as 
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23| 



S'S\A^)=I\Aa), 



(1) 



with the complex symmetric matrix S defined by 



£ = exp(i7rz>) U exp{infi) 



(2) 



where v and (x are diagonal matrices with diagonal elements z/^v and /iA and is a unitary 



matrix given right below. This equation has the following interpretation [2l|. I^a) is the 
set of angular coefficients in the "collision basis" (lA.ip of the wavefunction at the perigee 
of the trajectory. Eq. ([1]) means that at the quantized energy this wavefunction goes back 
onto itself when applying in order (from right to left) the following operations: 

1. a diagonal matrix with element exp(i7r/iA) = exp(— i7rA;/(47rL) A^) (see Eq. IA.4p . i.e. 
(half) a quadratic torque in the collision frame with an angle of tt, a strength of k/{A7TL) 
and an operator L'^Zq^ which brings the wavefunction from perigee to outgoing r = tq, 

2. a matrix Una which transforms from the "collision basis" (lA.ip labeled by A to the 
"free rotation basis" (\A.2\i labeled by A^, 



3. a diagonal matrix with element exp(i7rz/Ar) ~ exp(— ivr^A^), (plus a constant phase) to 



first order in A^ 



2l| . where and Tc are respectively the average periods of the electron 
orbit and of the free rotation of the core (see also the derivation below Eq. (|5])), i.e. 
approximately (half) a linear rotation with angle vr, a strength T^/Tc and an operator 
A^, which brings the wavefunction onto its apogee, 

and then applying the same operators in reverse order to bring the wavefunction back to its 
perigee. This means that at quantized energies the wavefunction is invariant (not merely 
unitary as at other energies) under the action of a matrix S which is a Quantum Poincare 
Map {21] in angular space from perigee to perigee. 

Notice then that by multiplying Eq. ([T]) to the left by exp(i7r/i) and reorganizing the 
result gives the equivalent equation 

exp(i27r/i) ?7*exp(i27rz>)f/|A^) = l\Al), (3) 

with 

A° = exp(i7r/i)A (4) 

which means with the same reasoning that the are the coefficients of the angular wave- 
function at the outgoing r = vq position in the collision basis, and that this wavefunction at 
quantized energy E is invariant by the quantum Poincare map given by Eq. ([3]). This is the 
equation we will employ to derive the kicked top model. 



4- Classical model: the stroboscopic map 

The preceding quantum theory has a classical counterpart [20| which is a succession of 
two different motions (see Fig. [11(b)): 

• When the electron is far from the core, their motions are not coupled, due to the 
rotational invariance of the Coulomb potential. The electron freely rotates around its 
angular momentum L fixed in space, and the core directed along M rotates around 
its angular momentum N perpendicular to it. When seen in the molecular reference 
frame, which OZq axis is M, the rotation of L is in retrograde direction. This is the 
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reason of the "anomalous commutation rules of momentum in molecular basis' 
The angle of rotation of this apparent motion at angular velocity Uc = d{BN'^)/dN = 
2BN during the time of an orbit of the Rydberg electron, which has angular velocity 
= d{-l/{2vl))/dvN = is 5<^L = -27r Te/Tc = -27r x 2BNv%. This is 

approximately only a uniform rotation of the average angle Sipi = — 27r x 2BJv^j when 
L ^ J, so that A^, which varies between J — L and J + L varies only slightly. 



During the collision the motions are coupled. L rotates around M by an angle 5ipL = 
—2nd^\/ OK = kA/L. Since the total angular momentum J = N + L (which are 
both well defined in classical mechanics during the collision) is conserved, this relative 
motion of L around M entails a "recoil" of the molecular reference frame whose precise 



value was computed in 
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B. The kicked top derived from the Rydberg molecule model 

1. Evolution operator and wavefunctions 

The preceding model is similar to the standard kicked top in that it displays a succession 
of quadratic torques around M separated by approximately pure rotations for L around the 
perpendicular axis N . But the non exactness of the pure rotation and the recoil motion 
(modifying N) during the collision step entails extra complexities (and extra interests arising 
from the explicit coupling of the angular momenta L and N). We now establish that for 
L <^ J, our Rydberg molecule model becomes a standard kicked top. Hence the kicked top 
can be considered to be a special limit of MQDT. 
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To this end we first compute (see Eq. (1A.9P ) for a generic value of as a function of 
the middle z/j, for N = J, which is related to total energy E hj E = BJ{J-\-l) — 1/ (2z/j). uj 
can be considered as another measure of total energy, in fact a linearized function of energy 
with ave rag e unity spacing in each series ( "unfolded" in the language of Random Matrix 
Theory 25|). 

= — 

v/l + 25z/2(Ar- J)(iV+ J+1) 



1 + 



Te N-J N+J+1 



Tc vj 



J 



-'- c 

with the average ratio of electron to core periods being given by 



Te _ UN=J _ dE/dJ 

T, We dE/duj 



2BJ uj. 



(5) 



(6) 



The limit in Eq. (E]) supposes that ^ and L (thus approximately N — J) remain constant, 
while uj oo. Thus according to Eq. ([6]) we must have 2BJ — )■ 0, which we suppose 
satisfied by 5 — )■ and J — )■ oo, thus also ~ J — t- oo, and (A^ + J + 1)/J — )■ 2. The 
important result is the linearity of with respect to J — N = Ml, which varies between 
—L and +L. Indeed in this limit is, to first order in L/J, equal to the projection of L 
onto J (obtained by developing A^^ = |J — Lp = J^ + L^ — 2L-J). J being invariant can 
be without loss of generality taken as laboratory OZ axis, so that Ml is the projection of L 
onto the laboratory OZ axis. 



Using t 
Edmonds 



le asymptotic formula for Clebsch Gordan coefficients (with the conventions of 
26l . A2.1]) when two angular momenta go to infinity, the remaining being finite. 



the elements of the transformation matrix become 



UmtK 



2:^-M,-A(0,7r/2,0) 



l)"+MM.(0,vr/2,0), 



(7) 



where the last expression makes use of 



271 ] or "passive" 
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, eqs. (4.2.5)-(4.2.6)]. Here 2) is the "standard" 
26[ rotation matrix, i.e. when rotating the reference frame with Euler 



angles a,/?, 7 while keeping fixed the quantum system: 



2)^ 



(Lm|e+^^-^^e+^'^^^e+^"^^ |Lm') 



(8) 



Explicit expressions are given in these textbooks 2m, |26|, l27|] , and we have written fast and 



accurate recursive programs programs for large momenta. Notice the + sign, and angle 



ordering opposites to the "active" point of view 28| (rotating the spin in a fixed reference 



frame), which is more common in kicked top works. The passive point of view is more 



"natural" in molecular works 



24|. 



Inserting the expressions ([5]) and ([7j) into Eq. ([3]) and moving the common factor 
exp(2i7rz/j) to the right hand side the quantization condition 



A/2 



^e-^^^2)^,^^,(0, -vr/2, 0)e^^-^*^-2)i^,jO, vr/2, 0)|i^) = e-'-^-^6^^,\Al) (9) 

A 

where we have defined the elements of \A\) by 

A" = (-1)-^-'^^°. (10) 

This change of sign takes into account that the molecular wavefunction contains a core 
part in addition to the Rydberg electron part (Eq. lA.ll) . We want an equation for the 
Rydberg electron only, and for L <^ J the core parts for different A differ only by this sign. 
To interpret the quantization condition, it can be noted that for a "passive" rotation, the 
transformation 

^k=E^A®AM.(0,7r/2,0) (11) 

A 

yields the coefficients of the outg oing wavefunction in a frame rotated around OY by 7r/2 
from the OZ axis (see e.g. 2^, Eq. (58.7)]), that is the OX axis. Alternatively notice 
that Eq. (I9l) is an eigen system of equations with eigenvalue exp(— 2i7rz/j) for a unitary 
map consisting from right to left: a frame rotation from OZ to OX, a pure rotation with 
parameter ^ along the new axis, back to original axis system (giving overall a free rotation 
along OX), and finally a quadratic torque with parameter k along the original OZ axis. 
Therefore this is exactly the evolution equation from kick to kick for the standard kicked top 



29| . Notice only the plus sign in the pure rotation term, which appears as the consequence 
of the apparent rotation in "wrong" sense when viewed in the molecular frame. Conversely 
one can say that this sign is correct in the laboratory frame, and that two other minus signs 
are a consequence of the map being described in the molecular rotating axis system. 
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III. ENTANGLEMENT IN A SINGLE KICKED TOP 



A. The kicked top as a compound system 

A kicked angular momentum J of the standard top can be considered as resulting from 



the composition of 2 J spin-1/2 subsystems so that J = Y^n=i^^ The spins s" must 

be in a state symmetric by permutation in order to generate the Hilbert subspace for the 
states \ JM) of the kicked top. Then the average for an individual spin s" projection on the 
axis i = x,y,z is identical for each n and related to the averages ( Jj) of the total angular 
momentum through 

(^D = (12) 

We will be interested in the entanglement between an individual spin s and the remaining 
subsystem containing the 2 J — 1 other spins. The average dynamics for s are obtained from 
the density matrix for the kicked top p{t) by taking the partial trace over the remaining 
spins. Note however that this step need not be done explicitly given that any spin-1/2 density 
matrix can be written as 

Ps{t) = \ + 2Y,{s,)s,, (13) 

i 

SO that the expansion coefficients are actually encoded in the averages. This is a considerable 
simplification relative to the modified kicked top. 



B. Entanglement generation and the quantum- classical correspondence 

1. Linear entropy 

We will quantify entanglement by computing the linear entropy 5*2 (t) associated with the 
reduced density matrix defined by 

S,{t) = l-TTpl{t). (14) 

For a pure state pi = ps and 5*2 vanishes, whereas for a maximally mixed qubit 5*2 = 1/2. 
By plugging Eqs. (fT2|) - (|T3l) into Eq. (fT4l) . we have 

S2{t) = \-^, + (Jy? + (Jf) ■ (15) 
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Hence entanglement depends on the averages of the kicked top angular momentum projec- 
tions; entanglement is maximal when all these averages vanish. 



2. Coherent states 



Angular momentum coherent states 30|] are the most suitable choice in order to inves- 



tigate the quantum-classical correspondence in the kicked top, and have consequently been 



employed from the early works onward 
angular momentum eigenstates by 



10|. These coherent states, given in terms of the 




|^^,0) = (l + tan^-)-^ > tan-e^<^ \JM) (16) 



are localized on the sphere, and in the present context they present the additional advantage 
of yielding an initial product state. Indeed, from the property 

{e,<p\j,\e,<p) = j„ (17) 

if p(t = 0) = 1^0) 0o) (^Oj 0o| "we then have {Ji(t = 0)) = Jj resulting in S2(t = 0) = 0. 

The most straightforward way of representing a coherent state and its ensuing evolution 
on the sphere is through the use of a Husimi distribution, which is precisely defined as the 
coherent state representation of the density matrix. For the standard kicked top in state {ip) 
this is simply given by 

h{e,<jy) = \{e,<jym\ (18) 

When lip) is itself a coherent state {tp) = |^^O)0o) the Husimi distribution is given by the 
overlap of the coherent states 

1(^,0 1^0, 0o)| =cos-'( ) (19) 

where x is the angle between the directions {0,(f)) and (6'o,0o)- This gives a distribution 
localized on (^O;0o) with an angular spread inversely proportional to J. 

In the rest of this work we will choose an initial coherent state centered on {6q, 0o) at 
t = and monitor the entanglement generation by computing 5*2 (t) that depends on the 
averages of the angular momentum projections through Eq. (|T5|) . 
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FIG. 2. Surface of section at "resonance" (see text). 
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FIG. 3. The linear entropy as a function of time when the initial state is a coherent state centered 
on (7r/2,0) (resonant case, k = 0.1). 

3. Classical distributions and averages 

Several aspects of the quantum-classical correspondence in the kicked top have been 



investigated 
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31| . Here we will only be interested in comparing the quantum averages 



(Jj(t)) to the corresponding classical averages Ji{t), when the initial quantum distribution is 
the coherent state 16*0, 0o) and its classical counterpart is a distribution of particles centered 
on (6'o,0o) and distributed on the sphere according to the right hand-side of Eq. ( fTOjl . For 
very short-time scales the classical and quantum averages are expected to be the same, 
though for longer times in typical cases at best a similar qualitative behavior can be obtained 
when the statistical averages are averaged out on relevant time scales in order to smooth out 
interference effects (such as quasiperiodic quantum revivals when the underlying classical 
dynamics is regular or the random superpositions when the underlying dynamics is chaotic). 

Our interest in the quantum-classical correspondence in the present context is to assess 
whether the configurations (kick strength, relative rotational and torsional frequencies, ini- 
tial state position) leading to efficient entanglement generation in the quantum top can be 
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FIG. 4. (Color online). Evolution of the quantum Husimi distribution (top row) and the analog 
classical distribution (bottom row). The plots show the short-time evolution of the distributions 
for times corresponding to the rise in S2{t) seen in Fig. [3] (resonant case, k = 0.1). The initial 
coherent state centered on the yellow dot in (a) and its classical counterpart in (c) is subjected to 
the torsional motion resulting in a spread of the distribution on a strip on both sides of the equator 
(the coloring reflects the intensity of the distribution) . 

related to different dynamical regimes of the classical top. For these purposes, a numerical 
comparison of the linear entropy [Eq. (ITB]) ] and of the corresponding classical expression 

C2{t) = \-^2{Jl + Jl + Jl) (20) 

is sufficient. We can further expect that the similar behavior of 5*2 and C2 is due to the 
similar behavior of the classical and quantum angular momentum averages. Contrarily to 
the situation of genuine 2-particle systems (in particular in our modified top), here C2 is 
totally unrelated to a measure of the non-separability of phase-space distributions that could 
classically play the role of 5*2 as a marker of classical mixtures Q]. Intuitively, it can be 
anticipated that starting from a localized state on the sphere, chaotic dynamics will tend 
to scatter the initial distribution over the entire sphere, leading to vanishing or small Jj for 
the three axes. This will lead to a maximization of C2{t), and should correspond to maximal 
entanglement generation. 
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IV. RESULTS 



We will compute numerical results for the quantum and classical tops with J = 10. This 
moderate value of the angular momentum is sufficiently low as to envisage the possibility 
of an experimental realization while already displaying the main features of the quantum- 
classical correspondence. We keep to the conventions introduced in Sec. II (i.e. free rotation 
around the OX axis and torsion along the OZ axis). 



A. Regularity at "resonance" 

In the Rydberg molecule, resonance refers to the electron period being an integer 



multiple of half the core period Tc (half is due to Kronig's symmetry of the core 33Q. This 



situation has observable consequences, appearing as clear zones in the spectrum [20|, and 
achieves high entanglement generation when the dynamics is regular ^\. In the standard 
kicked top limit, the free rotation becomes trivial: the classical dynamics is always regular 
(irrespective of k) and constrained to remain on the initial circle, while increasing k leads to 
an arbitrary separation between successive points on the circle. The corresponding surface 
of section is displayed in Fig. |2l 

Let us take k = 0.1 and an initial distribution centered in {6o,4>o) = (7i"/2,0) at the 
intersection of the x axis with the sphere (orange dot in Fig. [2]). The entanglement rate is 
shown in Fig. [3l We see that the linear entropy rises slowly and monotonically until it reaches 
its maximal value of 1/2. The classical quantity C2{t) follows strictly the same behavior. 
This behavior is due to the effect of the torsion on the distribution: in Fig. H] we have 
displayed the short-time evolution of the quantum (upper plots) and classical distributions. 
The initial coherent state, shown in Fig. Amoves to the left on the upper half of the sphere 
and to the right on the lower half of the sphere. This ensures that both {Jy{t)) and (Jzit)) 
keep their initial value of 0. {Jx{t)) on the other hand evolves from initial value {Jx{t = 0)) = 
J to {Jx(t)) = as the distribution stretches and encircles the sphere along the equator. 

This feature is readily understood by looking at the evolution of an initial classical distri- 
bution corresponding to the coherent state (bottom row of Fig.Hj): the particles far from the 
center of the distribution spread faster than those near the center and the first ones reach 
the opposite side of the sphere while the latter are still close to (6*0,00)- After a few kicks 
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FIG. 5. Top: Evolution of the linear entropy shown in Fig. [3] for longer times. Bottom: the 
analog quantity C2{t) for the corresponding classical system shows the same behavior, due to 
quantum-classical correspondence for the average of the projection J^. 

the distribution becomes approximately a uniform strip around the equator during which 
time Jx — and C2 — 1/2. For longer times the distribution relocalizes on the opposite side 
(7r/2, tt) (this is a purely classical effect) with almost equal to — J and spreads again. The 
corresponding behavior of 5*2 (t) for longer times is shown in Fig. |5l along with C2{t). 

Note that the time averaged entanglement rate before the first relocalization is extremely 
high, 5*2 ~ 0.46 (for longer times partial relocalizations proliferate and the average decreases 
to 5*2 ~ 0.43). It is interesting to compare with the case k = 10. The entanglement evolution 
is shown in Fig. El 5*2 (t) reaches the maximal value of 1/2 in only a couple of kicks, but 
periodically drops to significantly lower values (the time average is S2 ~ 0.43 identical to 
the k = 0.1 case). 

B. Mixed Phase-Space 

Let us now investigate the entanglement generation for the mixed phase situation Tg = 
0.95 and k = 5 (see the classical surfaces of section in Fig. |8]). We first take an initial state 
lying in the chaotic sea (Fig. [HI top panel). The evolution of 5*2 (t) is shown in Fig. [7](a). 
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FIG. 6. Top: S2{t) for the "resonant" case and A; = 10 (the initial state is the same coherent state 
as in Fig. [3]). 

The linear entropy reaches its maximal value of 1/2 after only a few kicks, but strong dips 
keep the time averaged value over the first 1000 kicks to 5*2 ~ 0.40. The evolution of 
the corresponding classical distribution is shown in Fig. |8] (a)-(c): the distribution quickly 
spreads over most of the chaotic sea, that nevertheless only covers a part of the available 
phase-space. 

Let us now take {6q, 0o) centered on a point lying in the large regular region as indicated 
in Fig. E] (lower panel). The linear entropy is shown in Fig. W[h). It immediately rises to 
5*2 = 0.45 and then oscillates wildly with a time average of 5*2 ~ 0.40. An inspection of 
the behavior of (Jj) shows that (Jy) and (J^) display quasi-periodicities typical of quantum 
revivals in regular systems whereas (Jx) follows on average the behavior of the classical 
average (see Fig. [9]). The evolution of the classical distribution, shown in Fig. |8]^d)-(e), 
indicates that most of the distribution stays within the island of regularity but spreads 
along the lines in the surface of section. The spreading is not perfectly uniform but the 
approximate symmetry of the distribution suffices to reduce significantly the averages Jj, 
which in turn increases C2{t) to values similar to the case when the classical distribution 
explores the chaotic sea. Note that in these two cases the time-averaged entanglement rate 
is exactly the same. 

C. Chaotic Phase-Space 

For Te = 0.90 and k = 15 there are no more islands of stability and the stroboscopic 
map is chaotic. The entanglement rate is shown in Fig. [TOl rising in a couple of kicks to the 
maximum bound. The time average of the linear entropy is 5*2 ~ 0.47. The quantum and 

17 



100 



200 



(b) 0.5 



c8'0.25 





100 200 

t (kicks) 



FIG. 7. S2{t) for the case Tg = 0.95, k = 5 when the initial coherent state is chosen to he in the 
chaotic sea [(a)] or in the central regular region [(b)], on the points shown on the surface of section 
in Fig. El 



classical averages of the angular momentum projections {Ji(t)) and Ji(t) oscillate around 
zero with a small but far from negligible amplitude (about ±2), thereby explaining the dips 
visible in the linear entropy. 



V. DISCUSSION AND CONCLUSION 

In our previous studies regarding the correspondence between entanglement generation 
and the underlying classical dynamics in Rydberg molecule top 5|, |9|], we had found that the 
linear entropy was correlated with the classical diffusion in regions of phase-space leading 
to inelastic scattering between the angular momentum of each of the two particles. Our 
conclusion was that chaotic classical dynamics tended to favor inelastic scattering, but was 
by no means a necessary condition, given that regular dynamics could achieve in certain con- 
ditions more efficient inelastic scattering (resulting in higher entanglement in the quantum 
system) . 

In the limiting case yielding the standard kicked top in which the single angular momen- 
tum is regarded as a collection of entangled spins, the correlation variable is considerably 
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FIG. 8. (a)-(c): Evolution of a classical distribution for the case Te = 0.95, k = 5 when the initial 
distribution is the analog of the coherent state centered on the green dot lying in the chaotic sea 
(whose linear entropy was plotted in Fig. [7]^a)). (a) shows the initial situation, (b) the evolution 
after 15 kicks, (c) the evolution after 75 kicks, (d)-(f): Same as above with the initial distribution 
being the analog of the coherent state centered on the yellow dot lying in the regular region (the 
linear entropy was plotted in Fig. El^b)). (d) shows the initial situation, (e) the evolution after 4 
kicks, (f) the evolution after 75 kicks (most of the evolution follows clockwise an elliptical motion 
and returns near the initial point). 

simpler: it only involves the sum of the angular momentum projection averages. The linear 
entropy given by Eq. ( [T5|) has a direct classical counterpart C2(t) given by Eq. (!20!) . Chaotic 
dynamics tends to maximize C2(t) as the initial delocalized distribution spreads all over the 
sphere. However we have seen that regular dynamics can also achieve high values of C2{t) by 
taking advantage of the symmetry of the evolving distributions. In the regular "resonant" 
case (Fig. ^ the time averaged entanglement before the first relocalization is equivalent to 
that of the entirely chaotic top of Fig. [101 actually contrarily to the chaotic case, in the 
resonant kicked top C2{t) sticks for a period of several hundred kicks to its maximal value 
of 0.5 without any dip. Hence the resonant top is an interesting candidate to control and 
achieve the highest degree of entanglement (the main drawback being the slow rise, which 
may conflict with decoherence scales in practical applications). 

The mixed phase-space top also shows interesting properties: the linear entropy when 
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FIG. 9. The averages of the angular momentum projections Ji{t) for a kicked top with Tg = 0.95 
and k = 5 are plotted for the quantum top (top row) and the corresponding classical top (bottom 
row) when the initial distribution is in the regular region. 



the initial state lies in the classically regular or chaotic regions shows a similar behavior, 
although the underlying classical dynamics is radically different, as portrayed in Fig. |8l The 
initial classical distribution in the chaotic sea mostly spreads within the available phase- 
space region, inducing a rapid rise in C2{t) (which is not maximal because a large portion 
of the sphere has regular features). For the regular initial distribution case the spread takes 
place essentially within the regular region. The averages shown in Fig. [9] indicate that Jy 
and Jz are not zero (as they would be for a uniform distribution) but oscillate instead with a 
small amplitude for and a larger one for Jy. Still, given that Jx itself oscillates around its 
small initial value, the resulting C2{t) is large, of the same order of magnitude as when the 
distribution was confined to the chaotic sea. Note that the oscillations of (Jy) and (J^) in 
Fig. ini that are typical of recoherences and revivals when the corresponding classical regime 
is regular cancel out when summed in the expression of 5*2 (t) and are therefore not visible in 
the entanglement rate. Here instead the time dependence of 5*2 (t) in Fig. [Tl^a) (initial state 
in the classically chaotic region) displays an aspect reminiscent of recoherences. 

The present findings do not disprove earlier results [xsj on chaos and entanglement 
in the kicked top in the sense that generically regular classical dynamics will tend to be 
correlated with lower quantum entanglement than when the classical dynamics is chaotic. 
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FIG. 10. 5*2 (t) for the case = 0.90, A; = 15 corresponding to a classically chaotic phase-space. 



Notwithstanding we have given explicit illustrations in the kicked top indicating that this 
generic behavior is not universal. The reason is that the linear entropy is correlated through 
the quantum-classical correspondence to the quantity C2{t) introduced above, which in turn 
depends on the classical averages. C2(t) does not depend simply on the global dynamical 
regime but on the details of the initial distribution and on its specific dynamical evolu- 
tion. In the examples we have given high values of C2{t) were obtained for distributions 
evolving through regular dynamics by appropriately choosing the localization of the initial 
distribution. Note that this feature can be enhanced by choosing multiply localized initial 
distributions (eg a sum of a couple of coherent states each centered on a different point on 
the sphere) so that the ensuing regular dynamics minimizes the averages entering C2{t) or 
5*2 (t). We therefore conclude that as was already seen on other systems 

QQ,™. general 

claims linking chaos and entanglement in the standard kicked top should be made with care, 
as they are not universally valid. 



To sum up, we have introduced the standard kicked top as the limiting case of our two- 
particle kicked top modeling Rydberg molecules employed in earlier works and investigated 
dynamical entanglement in the standard top as a function of the underlying classical dynam- 
ics. By linking the marker of the entanglement rate to a classical function depending on the 
angular momentum projection averages we have seen that the entanglement generation in 
the quantum kicked top depends on the specific details of the underlying classical dynamics 
rather than depending generically on the global properties of the classical regime. 
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Appendix: Rydberg molecule: Quantum solutions 



The solutions of the quantum problem are obtained by expressing the wavefunction in 
terms of two different angular bases in order to take into account the short-range phase- 
shifts (determined in the region near the core) and the asymptotic boundary conditions (for 
the radial variable going to infinity). The two expressions are matched at some radial value 
r = To for which both decompositions are valid {tq lying actually near the core). 

r < To Collision or molecular basis.: In this region the Rydberg electron is tightly 
bound to the core by a strong anisotropic potential, so that we have ordinjiry molec- 
ular diatomic eigenfunctions (Hund's case (b) in molecular nomenclature jsS]), where 
only the total angular momentum J = N + L is conserved (but not separately L and 
N) so that the good quantum numbers are its modulus J, its projection onto the 
OZ laboratory axis Jz = Mj, and its projection onto the axis of OZq = M of the 
core Jzq = Lzq = A. The last equations suppose that the core itself has no angular 
momentum along its own axis, i.e. is in a S state, so that the only angular momentum 
along the axis M is that of the Rydberg electron. To be definite, we choose the 



conventions of the original work of Fano 16|], namely 



The (quantum) molecular reference frame is obtained from the laboratory refer- 
ence frame by a rotation of Euler angles ipM, (with the "quantum" conven- 



tion that the second rotation is around the intermediate OY axis) 
• there is no extra overall phase factor 

so that the angular part of the basis wavefunctions is 



24 



26 



xi^iJ^(^^,</);,M) = Yi:{e',,<f>',)\l^^^i,,^{o,eM,^M) (a.i) 



where 6'^, ip'^ are the polar angles of the Rydberg electron in the molecular frame, and 
where the D are given in Eq. ([8]). 

r > ro Coulomb or laboratory or free rotation basis.: In this region the interaction 
potential is approximated by a pure Coulomb 1/r potential which is rotationally in- 
variant. The angular momenta of the core N and of the electron L are separately 
invariant in laboratory space, as is their sum the total angular momentum J. The 
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angular part of the basis wavefunctions is thus obtained by adding the two partial 
angular wavefunctions with a Clebsch-Gordan coefficient: 



AM^ 



An 

(A.2) 

the first index of 2)^ is because the core is in a S state. 

Quantization is obtained by demanding that the wavefunction go to zero at both ends 
r — )■ and r — t- oo. 

inner region r <= tq.: For a full fiedge solution, the radial equation is integrated outwards 
from r = to r = tq with the full anisotropic potential. This is done for each value of 
the total energy E and each value of the good quantum numbers, which in this range 
are according to Eq. (lA.ip L, J, A, Mj. In very serious cases one takes into account 



that for very short range the interaction is more complex than a mere potential. In 
any case at r = tq the radial equation reduces to a second order equation in the 
isotropic 1/r Coulomb potential, which local solutions are known. They are a linear 
combination of the regular and irregular solutions at origin r = of this problem as 

fAir) = (s(Ee, r) cos(7r/iA) + c(Ee, r) sin(7r/iA)) (A.3) 

where is the Rydberg electron energy {not the total energy E). It depends on only 
one parameter which depends itself on A, L, J (but not Mj). It depends also on 
Ee, but very slightly so for highly excited levels, because near r = tq the attractive 
Coulomb potential is much greater than the splitting of the rotational energies of the 



core, which are of the order of the sp 
tend to zero at the ionization limit 



itting of the higher energy electronic levels which 
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16| . Frequently, instead of computing /ia by 



radial integration, it is taken as a parameter to be adjusted to experiments. Due to 
Kronig's symmetry 33| by reflection on the molecular core axis, it must be a quadratic 



function of A [20|]. The simplest case is 

plus an unimportant constant, where is a classical parameter defined below. This 
is what we do always in this paper. 

23 



outer region r >= tq.: In this region up to oo the solution of the radial equation is of the 
same form as Eq. flA.Sp . but here the good quantum number is instead of A, so that 



it writes 

/jv(r) = {s{E„ r) Cn + c{E,, r) d^) (A.5) 
Furthermore, here Ee is known a.s E^ = E — BN{N + 1), where B is the rotational 



constant of the core, according to standard molecular conventions 



33|. 



matching at r = tq.: These two solutions are two developments of the same wavefunction 
on two different angular basis. Matching involves an overlap matrix 

^(LJ) _Y^^(LJ) friLJ) . . 



TV 



which is proportional to a Clebsch-Gordan coefficient 16| : 



f/Jf/) = {L - A, JA\Ljm){-iy-^+^, (A.7) 

Fano [16] gives also an extra symmetrization on ±A due to the symmetry of the 
potential by reflection onto the M axis. 

quantization.: It is obtained by demanding that the growing Coulomb solution at r — oo 
be zero, which is obtained if energy E is such that Eq. (lA.Sp reduces to 

/Ar(r) oc - cos(7rz/Ar) s{Ee, r) + sin(7rz/Ar) c{Ee, r), (A. 8) 

where un is the principal quantum number of the electron in channel A^ given in atomic 
units by 

E = BN{N + 1) + Ee = BN{N + 1) - ^ (A.9) 
The final quantization set of equations is 

Una sin(7r(z/^ + /iA))AA = ; N = J - L ■ ■ ■ J + L (A.IO) 

The standard way to solve it is to look to zeros of the preceding determinant (and thus non 
zero values of the coefficients A^) while varying total energy E, and thus un and possibly 
/iA which depend on E. 
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